Random Resonators and Prelocalized Modes in Disordered Dielectric Films. 
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Areal density of disorder-induced resonators with a high quality factor, Q > 1, in a film with 
fluctuating refraction index is calculated theoretically. We demonstrate that for a given kl > 1, 
where k is the light wave vector, and I is the transport mean free path, when on average the light 
propagation is diffusive, the likelihood for finding a random resonator increases dramatically with 
increasing the correlation radius of the disorder. Parameters of most probable resonators as functions 
^Lj ' of Q and kl are found. 
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PACS numbers: 42.25.Dd, 42.60. Da, 73.20.Fz 



Introduction. Recent discovery of the coherent lasing from various disordered materials adds a new dimension to 
the conventional physics of light propagation in multiply scattering media. In the experimental works it was 

demonstrated that, above a certain excitation level, the emission spectra of ZnO powders [0,^), conjugated polymer 
films |l and dye-infiltrated opals j| exhibit a sequence of extremely narrow peaks, their widths being limited by the 
spectrometer resolution. Until recently, this finite resolution left some room for doubt as to whether the observed peaks 
indicated true lasing Q. However, the latest experiments on photon statistics |6|,[?J have unambiguously established 
the coherence of the emitted light, thus proving conclusively that the underlying mechanism of random lasing in 
Refs. involves the amplitude (coherent) rather than power (incoherent) feedback. (The latter is known to occur 
via the diffusion process of the light intensity, as proposed long ago by Letokhov Q and observed first by Lawandy 
|)j and later in other numerous experiments.) The origin of coherent feedback, responsible for the lasing observed in 
Refs. 0-0], is a subject of controversy and debate. Such a feedback would naturally emerge if the light were localized. 
i^i - However, the coherent backscattering measurements Jl|,|l0|-p^| carried out in parallel with the analysis of the emission 
spectra rule out this possibility. Firstly, the values of kl extracted from these measurements turn out rather large. 
, Secondly, the onset of the Anderson localization manifests itself in the rounding of the top of backscattering cone 
JL5| p^t . No such rounding was observed in the experiments JUJloHhj ■ 
^*T* , In the early work of Cao et. al. (l],[l(],[ll|] , it was argued that, even in diffusive regime of light propagation, a 
photon, scattered off a certain sequence of impurities (grains of, roughly, lOOnm size) can move on a closed loop, 
which thereby serves as a laser resonator. The picture of such ring cavities put forward in Refs. [p] JIc| , p^|l was based on 
the experimental evidence |L8| that the recurrent scattering processes contribute to backscattering albedo. However, 
since in each scattering act most of the energy gets scattered out of the loop, an unrealistically high gain would be 
required to achieve the lasing threshold condition for such a loop. This point was particularly emphasized by van 
Soest 0, who also argued that "impurity loops" are likely to generate a broad frequency spectrum rather than 
isolated resonances. 

Certainly the picture of random cavities representing a certain spatial arrangement of isolated scatterers is too 
naive. This, however, does not rule out the entire concept of disorder-induced resonators. Although sparse, the 
disorder configurations that trap the light for long enough time can occur in a sample of a large enough size, and a 
single such configuration is already sufficient for lasing to occur. Therefore, under the condition kl 3> 1, which implies 
O ■ that overall scattering is weak, the conclusion about the relevance of random cavities can be drawn only upon the 
quantitative calculation of their likelihood. This is the subject of the present paper. 

We start with a remark that the issue of random resonators for light waves has its counterpart in the electronic 
transport. In particular, the modes with anomalously low losses are analogous to the so called prelocalized electronic 
states in diffusive conductors that are responsible for the long-time as ymp totics of the current relaxation. Theoretical 
study of these states was launched more than a decade ago (see Ref. [ ^9| ) and later renewed in Ref . |20) . The results 
obtained to date are summarized in the review pl[ . Following this analogy we dub the modes of random resonators 
with anomalously high quality factor, Q, as prelocalized modes. 

The principal outcome of our study is that for a given M > 1, the probability of formation of a high-Q random 
resonator depends crucially on the size of the scatterers, or, more precisely, on the correlation radius of the disorder, 
R c . Similarly to the treatment in Refs. 11 1 we restrict our consideration to the two-dimensional case (a disordered 



film). Regarding the geometry of a random resonator, we adopt the idea proposed by Karpov in Ref. J22| for trapping 
the acoustic waves in three dimensions. According to Ref. |p2| , the fluctuations responsible for trapping are the 
toroidal inclusions with reduced sound velocity. Correspondingly, in two dimensions, a random resonator represents a 



1 



ring-shape area (see Fig. la) within which the dielectric constant is enhanced by some small value e\ (compared to the 
background value eo). Then such a ring can be viewed as a waveguide that supports the modes of a whispering-gallery 
type. Due to the azimuthal symmetry, these modes are characterized by the radial index and the angular momentum, 
to. Denote by M m {kl, Q) the areal density of resonators with quality factor Q in the film with a transport mean free 
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path I. Here k — e ko, where fco = 27r/A, A stands for the wavelength in vacuum. Obviously, in the diffusive regime, 
kl > 1, the density M m {kl, Q) is exponentially small for Q> 1. In this domain Af m (kl, Q) can be presented as 

Nm{kl,Q)=N Q e- s ^ kl ^ , (1) 

where No is the prefactor. Now we can "quantify" the above statement about the dependence of the resonator 
likelihood on R c 

S m (k R c > 1) = <S>(e{ /2 k Q R c ) 
S m (k Rc « 1) 7rV2( e J/ 2 fc iJ c )3 ' 

The dimensionless single-parameter function $ is shown in Fig. 2. It is seen from Eq. (|J) that S m falls off rapidly 
with increasing R c . In the domain k R c > 1 but e\^ 2 k$R c < 1 we can set $ w 1, so that we have S m oc (koR c )~ 3 . 
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For larger R c , when koR c ^> 1, the function $ behaves as <I>(it) oc u. In this domain S m decreases slower with 
R c : S m oc {k$R c )~ 2 . We emphasize again that Eq. (||) applies for a given kl value, so that the decrease of S m with 
R c leaves the backs cattering cone unchanged. Below we will demonstrate with rigorous calculations that the value of 
S m (k R c <C 1) is really large 
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S m (k R c « 1) = 2 ^—j kl\ InQI , (3) 

which makes formation of random resonators with Q > 1 practically impossible for point-like scatterers. On the 
contrary, for k$R c > 2 ring-like fluctuations of the dielectric constant (see Fig. lb) make realistic candidates for 
high-Q random resonators. 

Derivation of Eq. (j^j. Denote by P the probability of fluctuation of dielectric constant, Se(r). Assuming the 
fluctuations to be Gaussian with r.m.s. A, we have 

InP = / 1 drdr 1 8e{r)5e{r')n{r - r'), (4) 

where the kernel n(r) is related in the usual way J dr'n(r — r')K(r' — ri) — S(r — r^), to the correlation function, 
K(r), defined as 

(Se(r)Se(r')) = A 2 K(r-r') . (5) 

We search for fluctuations of a type shown in Fig. 1, i.e., Se is azimuthally symmetric (depends only on the radius, 
p) and is non-zero within the relatively narrow ring of width w <ti po, where po is the radius corresponding to the 
middle of the ring (see Fig. la). For such a fluctuation the wave equation for field distribution corresponding to the 
angular momentum, to, reduces to 

d 2 Xrn 

Lx m = -pr + Se k 2 x m = £i k 2 x m , (6) 

where x = p — po- The "eigenvalue" in the r.h.s. of Eq. (||) is defined as ei = (^j^^j ~ eo- It is independent of x by 

virtue of the condition po 3> w. The next step is to find the "most probable" distribution 8e(x) for a given t\. This is 
done within the standard optimal fluctuation approach [ p3|j2^ ] , which prescribes minimization of auxiliary functional 
^{^e} = InP — X(x m Lxm), where A is the Lagrange multiplier. This yields 



Se(x) = J dxiK (x - xx)x m {x\) , (7) 

where we have set A = 27r ( oo/A 2 fc ! . The above expression reminds the standard result of optimal fluctua- 
tion approach for a correlated random potential |23} ]. The only difference is that, due to the angular integra- 
tion in Eq. (||), the kernel in Eq. ([?]) is given by a function K Q , which is related to the correlator, K{r), as 
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Kq(xi — X2) = dyK(yJ (xi — X2) 2 + y 2 )- A natural x-scale for the eigenmode \ m is (ej^fco) 1 . Thus, we 
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present the outcome of the optimal fluctuation approach in terms of a dimensionless variable 

Sm = e ^ k 2° A 2 J I dz l dz 2Xli(zi)xin( z 2) K o(zi - Z 2 ) . (8) 

The dimensionless equation for the function Xm reads 

d 2 Xm{z) , Xm(z) 



dz 2 ij 2 h 



dziK Q (z - zi)x m (zi) = Xm(z) ■ (9) 



To proceed further we have to specify the correlator, K. We have chosen the Gaussian form K (p) = exp f— p 2 /i? 2 ). 
The form of the function x-m = A(R C ) exp[— 7(i? c )z 2 ] allows to cover the entire range of correlation radii, from "white 
noise" (R c — > 0) to the limit of a smooth disorder (koR c 1). Indeed, for large R c this form becomes exact. In 
the opposite limit, R c — * 0, using the above trial function instead of exact solution Xm °c l/cosh(z) leads to the 
overestimate of S m by a factor (jr/3) 1 / 2 w 1.023. The parameter A(R C ) and 7(i? c ) of the trial function can be found 
analytically. Evaluation of Eq. (g) reduces to the Gaussian integration which for a given m « e 1 / 2 fcoPo yields 

The analytical expression for the function introduced above in Eq. (||), and shown in Fig. 2, is the following 

^)^(Itct) 2 [» 2+ > + ^»] ,/2 - 

Recall that we are interested in the density of random resonators at a given value of kl. The remaining task is to 
express the transport mean free path in terms of A and R c . With K(p) — exp (— p 2 /i? 2 ) the expression simplifies in 
two limits 

kl = n t.^2 for k o R c « 1 ; ( 12a ) 

kl = 4 ^J a ^° for k R c »l . (12b) 



Note, that Eqs. ( 12a ) , ( |l2b| ) can be cast into the conventi ona l form I cx 1/ na, where n is the concentration of scatterers 
and a is the transport cross section. In particular, Eq. (|!2aD follows fr om t he two -dimensional version (er cx R^k^) of 



the Rayleigh scattering formula for a. Combining Eq. (|10|) and Eqs. (12a), (12b), we arrive at Eq. (|2|). In the limit 
el^k^Rc ^ 1 Eq. ( |Io| ) has a simple interpretation. Namely, in this limit S m can be rewritten as S m ~ (ei / A) 2 (po / R c ) . 
The first factor comes from the Gaussian probability to have a fluctuation t\ within a square of an area ~ R 2 . The 
second factor accounts for the number of squares needed to cover the ring area A m ^ p$R c . 

l li 

Q-factor. It is obvious that at large distances, p po, the behavior of Xm is oscillatory, Xm oc exp(ie kop), 
manifestingthat the waveguided mode of a ring, being prelocalized, has a finite lifetime. In other words, the eigenvalue, 
ei, in Eq. (||) has an imaginary part, ei, due to evanescent leakage. The quality factor is inversely proportional to ?i, 
namely, Q = eo/ei- The leading contribution to e\ comes from the region of a width, d, adjacent to the waveguide 
(see Fig. la). To find d we have to take into account that the r.h.s in Eq. (||) does in fact depend on x. This is because 

the precise form of the r.h.s is not (^^j ~ e oko, but rather ( p "^ x — eofcp. In the region outside the waveguide, 

when x ^> l/e^fco (but still x <§; po) Eq. (||) takes the form 



eifco (l - | ) Xm , (13) 



dx 2 

where the width of the decay region is given by d — eik 2 pf ) /2m 2 . Equation ( p"3"| ) is of Airy- type. Semiclassical 
calculations with exponential accuracy yields for the rate of evanescent leakage t\ oc exp[— (2m/3)(ei/eo) 3 ^ 2 ], and 
hence 



3 



o / \ 3/2 



Substituting ei from Eq. (14) and A from Eq. ( 12a ) into Eq. (|10|) we arrive at Eq. (|3j). The above derivation relied 
on the assumptions po 3> d and d 3> w. These assumptions are justified within the following domain of the quality 

factors m ^> InQ ^> max jl, eQ^ 4 (/c -R c ) 3 / 2 m _1 / 2 |. 

Discussion. Equation (0) quantifies the effectiveness of trapping of light in a random medium with point-like 
scatterers. It follows from Eq. @ that the likelihood of high-Q cavity is really small. Indeed, even for rather strong 
disorder, kl — 5, the exponent, S m , in the probability of having a cavity with a quality factor Q = 50 is close to 

S m = 120. We emphasize that in two dimensional case under consideration, this exponent does not depend on m and, 

l li 

thus, on the cavity radius po = m / e o ^o- To estimate the degree to which finite size of sca ttere rs (~ R c ) improves 



the situation we choose k$R c ~ 2, which already corresponds to the limit k$R c 3> 1 in Eq. ( |12b| ), but still allows to 
set $ = 1. Then for Q = 50, kl — 5 we obtain S m ~ 1.1, suggesting that the resonators with this Q are quite frequent. 
In the latter estimate we have set eo = 4. 

A natural question to address is how large a value of Q can be achieved for a given kl and k$R c . To address this 
question we inspect the argument u = e\^ 2 koR c of the function $. With the use of Eq. ( fi~4| ) it can be presented 



in the form u = e^ 2 koR c 



31nQ 
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. Since $(u) increases monotonically (see Fig. 2), the latter expression suggests 



Ira 

that Q can be increased at the expense of larger m-values. In the example considered above, in order to keep $ 
smaller than, say 1.5, m should be bigger than 50. However, due to slow dependence u oc m -1 / 3 , we get rather 
small value S m rs 3 for m as small as m = 15. Certainly, the allowed values of m are limited from above. This 
limitation originates from "vulnerability" of waveguiding to the fluctuations of the dielectric constant around optimal 
ring-like distribution. The dangerous fluctuations are those that enhance the evanescent leakage. Note, that these 
fluctuations do not affect the main exponent S m in the density of resonators. It is obvious that the bigger is the area 
A m = 2npod oc m 4 / 3 , responsible for evanescent leakage (see Fig. la), the harder it is to "protect" the waveguiding. 
Since the fluctuations 6e(r) have the spatial scale R c , the probability that the waveguiding "survives" can be roughly 
estimated as exp(— A m /R^). The condition that the exponent A m jR\ does not dominate over the principal exponent 
S m can be rewritten as m < (e /ci) 3//5 ( m Q) 2 ^- F° r the example kl = 5 and Q = 50, addressed above, we get m < 10. 
Rigorous calculation of the "survival probability" is the problem of the same complexity as calculation of the prefactor 
in the functional integral |25| . 

Conclusion. In the present paper we provided a quantitative theory of random resonators that substantiates the 
intuitive image Jl],[l0| of a resonant cavity as a closed-loop trajectory of a light wave bouncing between the point-like 
scatterers. The intuitive picture in |l],[wj assumed that light can propagate along a loop of scatterers by simply being 
scattered from one scatterer to another. Such a picture, however, is unrealistic due to the scattering out of the loop 
0. We have demonstrated that the scenario of light traveling along closed loops can be remedied. In our picture the 
"loops", i.e. the random resonators, can be envisaged as rings with dielectric constant larger than the average value. 
On a microscopic level these resonators correspond to certain arrangements of scatterers (grains). The main point, 
however, is that a resonator acts as a single entity: only the coherent multiple scattering of light by all the scatterers 
in the resonator can provide trapping. We also point out that correlations in the fluctuating part of the dielectric 
constant (due to finite grain size) highly facilitate trapping. 

The effectiveness of light trapping is expressed by Eq. (j|). This expression describes the statistics of the quality 
factors which determines the distribution of the threshold gain for random lasing. Our consideration pertains to the 
passive disordered films, in the sense, that we neglect the effect of gain on the spatial distribution of the light intensity 
P, pl]j2^j27[ | . Random resonators considered in the present paper are sparse, so that there is no spatial overlap between 
the modes of different resonators, the situation opposite to that considered in Refs. p8p9| ]. We have also treated 
scatterers as frequency-independent fluctuations, 6e(r), of the dielectric constant. The entirely different scenario of 
the collective mode formation emerges for resonant scatterers |p0[| . 
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FIG. 1. (a) The structure of a two-dimensional resonator is illustrated schematically; only a half of the ring-shaped waveguide 
(blank region) is shown, (b) Optimal fluctuation of the dielectric constant, 5e(p) (solid line), and the corresponding field 
distribution (dotted line) are shown. Dashed line outside the shaded region of a width, d, illustrates the evanescent leakage. 
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